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The renormalization-group (RG) method in the theory 
of turbulence is based on the stochastic Navier-Stokes 
equation with a Gaussian random force Q, |^ |^ . One of 
the central problems in this approach is calculation of the 
Kolmogorov constant C, i.e. the dimensionless amplitude 
in the scaling law (Jj 

S2{r) = C(fr)2/3 (1) 

expressing the dependence of the second-order structure 
function 5*2 (f") on the relative distance r in the inertial 
range r^ ^ r L. Here, L is the external length of 
turbulence, the dissipative length and £ the energy 
injection rate per unit mass (which, in the steady state, 
coincides with the dissipation rate). 

Several attempts have been made in the past to solve 
this problem - 0|, but they all suffer from ambi- 
guities in connecting model parameters and observable 
quantities. As a consequence, there are significant dis- 
crepancies in the predicted numerical values for the Kol- 
mogorov constant (the spread is about a factor of two). 
In this Letter we analyse reasons of this unsatisfactory 
situation and present results of a calculation based both 
on an expression of C in terms of universal quantities 
and account of additional singularities arising near and 
below two dimensions. Rather unexpectedly, the analysis 
reveals that these singularities have a major effect on the 
numerical values of observable quantities well above two 
dimensions. We also show that a partial summation of 
these singularities is possible and significantly improves 
the numerical value obtained for C. To assess properties 
of the expansion produced within the RG approach, we 
have carried out the calculation in the two-loop approx- 
imation (the results of Refs. - 0| were obtained in 
the one-loop approximation). 

In the RG approach to c?-dimensional turbulence a 
powerlike correlation function of the random force is of- 
ten used: (/ /) ~ Dq /fc4-d-2e = df{k). In the RG frame- 
work various quantities may be calculated in the form 
of an e expansion which subsequently must be extrapo- 
lated to the physical value e = 2. For some important 



quantities the e expansion breaks off, which for the func- 
tion S2(r) yields the Kolmogorov exponent 2/3 [as in 
Eq. iCJ] at e = 2. To find the Kolmo gorov constant the 
amplitude of this function has to be calculated, which, 
however, can be done only approximately, because its e 
expansion does not break off. In calculation of the ampli- 
tude, apart from technical difficulties at two-loop order, 
a principal problem arises as well: the answer for S2{r) 
has to be expressed in terms of the energy injection rate 
£ [as in Eq. (^] instead of the parameter Dq of the pow- 
erlike forcing function. This problem has been treated in 
different ways in Refs. - which has led to different 
one-loop values of C. 

In Ref. (see also the connection between Dq 

and £ was sought with the aid of the exact relation 

In the unphysical region e < 2 this integral has to be 
cut off at wave numbers of the order of A = r^^. At 
fixed £ this procedure introduces, first, dependence on 
£ of the form Dq ^ (2 — e) in Dq (which has to taken 
into account in the construction of the e expansion) , and 
second, an ambiguity connected with the possibility to 
replace the upper limit A by cA with an arbitrary coef- 
ficient c. The first feature is rather natural, because the 
powerlike forcing df{k) ~ (2 — e) reproduces in 

the limit e ^ 2 the realistic forcing by infinitely large 
eddies: df{k) ^ <5(k). The second feature, however, in- 
troduces arbitrariness in the sought connection between 
Dq and £ through the coefficient c^^""*, which in turn 
renders the e expansion of _Do-dependent quantities am- 
biguous (in Ref. the simplest choice c = 1 was used ) . 
This is a reflection of the fact that the physical content 
of the theory remains unaltered when Dq is multiplied 
by an arbitrary function F{e) with F{2) — 1. 

Another way to fix the connection between Dq and 
£ has been used in Refs. Q - It amounts to the 

use of an exact relation (for the physical value e = 2 of 
the falloff exponent) which allows to connect the spectral 
energy flux with an integral of a third-order correlation 
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function, the latter being subsequently calculated in the 
form of an e expansion. The use of this relation in the 
unphysical region e < 2 is tantamount to a certain choice 
of the function F{e) mentioned above. 

Thus, the e expansion of the Kolmogorov constant in 
the model with the power like forcing is not unambiguous. 
Therefore, a better or worse agreement with the experi- 
mental value of C at one-loop level does not bear much 
meaning until a procedure for subsequent approximations 
has been pointed out and the stability of obtained results 
checked. On the other hand, since the real value of the 
expansion parameter e = 2 is not small, it is difhcult to 
expect good quantitative results without estimating - at 
least approximately - higher orders of the e expansion. 

In the model at hand, only quantities independent of 
Dq have rigorously unambiguous dependence on e (we 
will call them universal). Such quantities are, e.g., criti- 
cal exponents and dimensionless ratios of structure func- 
tions Sn{r), the skewness factor S = 83/82^^ in partic- 
ular. Calculation of universal quantities with the use of 
the RG method and the e expansion yields unambiguous 
results and cannot lead to such "paradoxes" as different 
one-loop values for the Kolmogorov constant. 

In view of this we have pursued the goal of find- 
ing a suitable universal quantity the physical value of 
which would be simply connected with the Kolmogorov 
constant, and calculating this quantity with the aid of 
the RG. The skewness factor S, connected with the 
Kolmogorov constant through an exact relation C = 
(—4/55)^/'^ might serve as such a quantity. How- 

ever, in the unphysical region e < 2 the structure function 
82 (r) in the model with the powerlike forcing correlation 
df{k) ^ k*^'^^'^'^ contains ~ at £ < 3/2 - an independent 
of r UV-divergent additive term ~ A^~^^^^ [for 83{r) 
this problem is absent, see below]. As a consequence, 
a straightforward generalization of the skewness factor 
S = 83/82^'^ to the region e < 3/2 becomes pointless, 
because the powers of r in this definition do not cancel 
due to the constant term in 5*2 (r). Therefore, as the de- 
sired universal quantity we chose the "nearest relative" 
of the skewness factor, the quantity 



_ rd82{r)/dr rd82(r) / dr 

" |^3(r)|2/3 = (-^3(0)2/3 



(3) 



independent of r in the whole region < £ < 2 which 
allows to find the physical values of S and C through the 
simple relations 



S = - 



-1 3/2 



3g(2) 



c = 



12 



d(d + 2) 



2/3 



Q(2). (4) 



The use of S'3(r) in constructing universal quantities is 
advantageous because it can be found exactly from the 
spectral energy balance for all £ < 2: 



3(d~ l)r(2-£) {r/2f'-'^DQ 
(47r)'i/2 r(d/2 + £) 



This expression allows, on one hand, to avoid calculation 
of graphs in construction of the e expansion for 83{r), 
and, on the other, confirms that passing to the physical 
limit £ ^ 2, in which r(2 — £) ~ 1/(2 — £), requires 
the dependence Do ~ a(2 — e) to arrive at a finite value 
of 5*3 (r) . The choice of the coefficient a consistent with 
ields the "4/5 law" of Kolmogorov: S'3(r) = ^^S-r 
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In the usual e expansion at d > 2 the universal quantity 
Q{e) © has the form 



Q(£)=£l/3Vgfc(d)£'= 



fc=0 



(5) 



The RG method allows to find successively the coeffi- 
cients of Qk{d) as a result of calculation of renormal- 
ization constants and scaling functions in perturbation 
theory (loop expansion). In Refs. - 01 only the one- 
loop approximation was used in the calculation of the 
Kolmogorov constant. The results of a two- loop calcula- 
tion with the aid of the relations Q , (|SJ have been quoted 
m Ref. fiej. For the one-loop contribution to Qo{d) in 
Eq. |3l an analytic expression for all d may be obtained: 



Qo(rf) = (l/3)[4(d + 2)]i/^ 



(6) 



The two- loop contribution Qi{d) gives rise to integrals 
which may be evaluated numerically for any particular 
values d. For d = 3 in the calculation of the Kolmogorov 
constant according to Eq. (El) the values C^^) = 1.47 
(one-loop approximation) and C^^^ — 3.02 (two-loop ap- 
proximation) were obtained. Although the two-loop cor- 
rection is not small, the recommended experimental value 
of the Kolmogorov constant C w 2.0 |j, |l7| turned out 
to be in between the values given by the two approxi- 
mations. Hardly any more could be expected in view of 
the fact that the value of the expansion parameter is not 
small. In what follows, we will show that the agreement 
with the experiment may be significantly improved by 
an approximate account of the high-order terms of the 
expansion ((SJ. 

Analysis of the dependence of the functions Qk{d) 
on the space dimension d shows that they have singu- 
larities at d < 2. In particular, Qk[d) ^ A"*"' for 
2 A = d — 2 0. This means that in the course of 
d tending to 2 the expansion ISJ necessarily will become 
" spoiled" , because the relative contribution of the high- 
order terms will grow without limit. In the present two- 
loop approximation this feature shows in that the ratio 
Qi{d)/Qo{d) in the limit d —>■ 00 (far away from all sin- 
gularities) is about 1/20 and monotonically grows with 
decreasing d assuming at c? = 3 the value ~ 1/2 of which 
the major part is brought about by graphs singular in 
the limit 2A = d — 2 ^ 0. This gives rise to hope that 
summation of leading A singularities in Eq. ISJ allows 
to improve quantitative results of the RG theory. 
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In the theory of turbulence the space dimension d = 2 
is exceptional from both the physical point of view (addi- 
tional conservation laws, inverse energy cascade) and the 
formal procedure of UV renormalization, because in the 
limit d — > 2 new divergences appear in the graphs of the 
perturbation theory. These divergences show in the form 
of poles in A in the coefficients Qn{d) for n > 1 in Eq. 
(0) . A consinstent procedure to remove these divergences 
by an additional renormalization has been developed and 
gives rise to a two-parameter e, A expansion In this 
Letter it is not possible to dwell on details of calculations. 
However, we want to point out the following principal is- 
sue. 

The use of the e, A expansion in the theory of turbu- 
lence was proposed in Ref. [l^ i whose author points out 
that an additional renormalization of the random force is 
required. This was carried out, following Ref. Q], by mul- 
tiplicative renormalization of the random force. There is, 
however, a major difference between the models of Refs. 
Q and |0| . In Ref. a local correlation function of the 
random force is considered (~ - model A, or ^ const 
- model B), wheras in Ref. the correlation function 
is nonlocal ^ k'^~'^~'^'^. From general theory of renor- 
malization it is known that counter terms may only be 
local (see, e.g. 113,01), which means that the renormal- 
ization adopted in Ref. is inconsistent. A consis- 
tent pro cedure, which we have used, was put forward in 
Ref. [131 • Our two- loop calculation allowed to confirm 
directly the general statement and show that it is not 
possible to renormalize the present model by multiplica- 
tive renormalization of the random force. We emphasize 
that, although this fact begins to show in the two- loop 
approximation, it renders the corresponding one-loop re- 
sult incorrect as well, since the use of the RG approach 
is based on the existence of certain relations between all 
terms of the perturbation expansion which break down 
in an incorrect renormalization. A detailed discussion of 
this issue we defer to a separate publication. 

In Refs. [1^ this two-parameter renormalization 
procedure (in one-loop approximation) was considered 
an alternative to the usual e expansion. We exploit it in 
a different manner - as a way to improve the expansion 
(|5|) by carrying out an approximate summation of the 
high-order contributions. 

To single out the leading poles, we express the coeffi- 
cients Qk{d) in the form 



(7) 



with a regular function 



(7fc(A)=^<Z«A'. (8) 

1=0 

Substitution of the expressions from Eqs.([7|l and (|SJ| in 



Eq. (O leads for the quantity Q to the representation 



Qie) 



a/3 



oc oc 



(9) 



k=0 1=0 

The £, A expansion corresponds to the asymptotic regime 
£ A ^ 0, A/e = const. Hence, the quantities {e/ A)^ 
in Eq. (0 are not considered small and the powers A' 
play the role of a formal small parameter. The quantity 
Q from Eq. ^ in the nth-order approximation (n > 1) 
is the series 



(10) 

fe=o 1=0 

which corresponds to an approximate calculation of the 
coefficients (fTj) of the e expansion jSJ with the account 
of n terms in the sum (|H| . For a RG calculation of the 
quantity qI."^ in the £, A-expansion scheme an n- 
loop approximation would be needed. 

Let us assume for the moment that we have carried 
out an n-loop calculation in the usual e expansion thus 
determining the following partial sum of the series © 

n-l 



i/3VQ,(d)e'=^Q, 



(n) 



k=0 



(11) 



and an n-loop calculation in the e, A-expansion scheme 
as well, hence having determined the quantity ^ of 
Eq. (|10|l . Then we have the possibility to amend the 
sum (|11() by an approximate contribution of all higher 
powers of not taken into account in Eq. (|ll(l . The 
required information of this contribution is contained in 

(n) 

the quantity Q^^. To obtain the improved value of Q 
we add the expressions ()10|l and (|11() . then subtract once 
the sum 



6Q 



(n) ^ ^1/3 



fc=0 1=0 

which enters twice in the sum of Eqs. (|10|l and Hll(l . 
Thus, we arrive at the following n-loop approximation 

+ Q^l - ^Q'"^ (12) 

for Q. Our two- loop calculation yields the result 



2(e + A)2£ 



Q 



e,A 



for the quantities , Q^^a "with the subsequent ex- 

pressions for 5Q^^^ , SQ'^^^: 
2 



3(2e-|-3A)2 
0.5181e 

-,(2) 



1/3 



-A 



(13) 



3 

^(1 



(28) 



1/3 ^ 

2e 
9A 



0.518l£ 



-A 



(14) 
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TABLE I: One and two-loop values of the Kolmogorov con- 
stant in the usual e expansion (Cg) and the double e, A ex- 
pansion (Ce.a); the contribution Cs in Eq. Q from the cor- 
rection (5Q'"' in Eq. I12L and the value Cejj from Eqs. 



n 






Cs 




1 


1.47 


1.68 


1.37 


1.79 


2 


3.02 


3.57 


4.22 


2.37 
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Calculating at = 3 the quantity Qe from Hll|) with 
the aid of ® and the value ^2(8) ~ 0.4748 found in 
and substituting the result together with Eqs. I|13|l and 
(|14|l in Eq. (|12(l we find the quantity Qeff in first and 
second approximation: Q^^j = 1-38, Q^^ff = 1-84. Sub- 
stitution of these values in Eq. Q at c? = 3 yields for the 
Kolmogorov constant and the skewness factor the values 
^eff = 1-79, Cff\ = 2.37, S^^fj = -0.33, = -0.22. 

In the Tablem [we have quoted for comparison the values 
of the Kolmogorov constant calculated according to Eq. 

in the first and second order of the usual e expansion 
(Ce), the double e, A expansion {C^.a), the contribution 
Cs in Eq. Q from the correction (JQ^"' in Eq. ifT^ 
and the value Ce// obtained from Eqs. I@J and (|12|l . 
In all the cases quoted the recommended experimental 
value of the Kolmogorov constant C ~ 2.0 lies between 
the values of the first and second approximation. How- 
ever, the difference between successive approximations is 
rather significant both in the e expansion and in the e, A 
expansion, let alone the leading terms of the e expansion 
of the latter. For the improved e expansion, i.e. for the 
quantity Cg// = + C^^a — Cs calculated according to 
Eqs. (|12|l and however, this difference is about three 
times smaller leading to a far better agreement with the 
experimental data. 

In conclusion, we have shown that a proper account 
of the "nearest singularity" in the coefhcients of the e 
expansion ij^l leads to a significant improvement of the 
results of the two-loop RG calculation at d = 3. We have 
analysed the effect ot this procedure at other d as well. 
It turned out to reduce significantly the relative contri- 
bution of the two-loop correction in the whole range con- 
sidered 00 > d > 2.5. At the same time this contribution 
remained large at d = 2, which we think to be an effect 
of singularities at the next exceptional dimension d = 1 . 

Obviously, the proposed procedure of approximate 
summation of the e expansion is applicable not only to 
the calculation of (3(e), but all universal quantities such 
as dimensions of composite operators. 
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